Entanglement of a qubit coupled to a resonator in the adiabatic regime 
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We discuss the ground state entanglement of a bi-partite system, composed by a qubit strongly 
interacting with an oscillator mode, as a function of the coupling strenght, the transition frequency 
and the level asymmetry of the qubit. This is done in the adiabatic regime in which the time evolution 
of the qubit is much faster than the oscillator one. Within the adiabatic approximation, we obtain 
a complete characterization of the ground state properties of the system and of its entanglement 
content. 
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I. INTRODUCTION AND DESCRIPTION OF 
THE MODEL 

The spin-boson model has been widely used to investi- 
gate the interaction between a two-level system, a qubit, 
and an harmonic oscillator environment, describing fluc- 
tuations of either electromagnetic or elastic origin 
The coupling of the qubit with each environmental mode 
gives rise to a progressive entanglement, leading to the 
decoherence of the qubit itself. 

This model has been largely employed in the weak cou- 
pling limit to explain noise effects in solid state devices 
which could be useful for quantum information process- 
ing 0|- It has been also applied to describe the coupling 
of such devices to quantum detectors. 

In the latter perspective, the strong coupling to a single 
bosonic mode has been analyzed in Ref. Q ■ This kind of 
"restriction" to a single-mode environment appears to be 
useful for the decoherence problem too, as recent exper- 
imental and theoretical works have attributed a promi- 
nent role to the coupling of superconducting Josephson 
qubits with spurious micro-resonator resulting from the 
presence of switching charged impurities residing in the 
tunnel barrier H 

In this paper, we analyze the case of a qubit strongly 
coupled to a slow resonator, working in the adiabatic 
regime. Our model is a generalization of the one em- 
ployed in Ref. Q, which turned out to also describe 
molecular Jahn- Teller effect 0. It can be also used to 
describe the coupling of a Josephson charge qubit to an 
electromagnetic resonator or to another (large) junc- 
tion working in the harmonic regime^, in the case of 
strong and off-resonant interaction. As it occurs for many 
solid state implementations, we assume, here, that the 
coupling can become so strong that the usual rotating 
wave approximation cannot be employed. 

Our aim is to characterize the ground state of the sys- 
tem and, in particular, to evaluate the amount of quan- 
tum correlation present in (that is, the "entanglement 



content" of) the fundamental level. If the presence of 
the oscillator is spurious and un-wanted, this "residual" 
entanglement can produce errors in the information pro- 
cessing performed by the qubit. An investigation of the 
entanglement in the case of a two state system coupled 
to an ohmic environment has been performed by Costi 
and McKenzie 0, by exploiting the equivalence to the 
anisotropic Kondo model. They were able to show that 
the entanglement entropy, for level asymmetry different 
from zero, reaches a maximum at smaller values of di- 
mensionless dissipation strength. 

In fact, the calculation of ground state entanglement 
has been used to characterize complex quantum many 
body systems, with particular emphasis on its connection 
to quantum phase transitions [iflllBllJillS- In our case, 
the system is bi-partite, and therefore there is no collec- 
tive behavior to be examined; but nevertheless, a kind of 
criticality is present, as in the massive limit for the oscil- 
lator (and for qubit working at degeneracy), two regions 
in parameter space exist, exhibiting completely separa- 
ble and entangled ground state, respectively Q- Further- 
more, a sharp increase from zero is found at the onset of 
entanglement. This has been interpreted as a quantum 
reminiscence of the bifurcation of the fixed point of the 
oscillator in the corresponding classical model 0. 

We show below that, within the adiabatic approxima- 
tion scheme, this behavior can be obtained analytically 
together with the leading corrections for a finite tun- 
nelling amplitude of the qubit. Our approach, however, 
is not limited to this region and we show that it can be 
used to systematically investigate ground state proper- 
ties and entanglement in a broader parameter range, as 
we can account also for the effect of level asymmetry. 

A qubit interacting with a single harmonic oscillator 
mode can be described by the Hamiltonian (in unit such 
that h = c = 1) 



A 



= («■'■ + a) 



■ ioc^ a 



(1) 
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where A is the transition frequency of the qubit, e is the 
level asymmetry, lo is the frequency of the oscillator and 
A is the coupling strength. 
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Hines et al. [g start their description from the case of 
frozen qubit (i.e., A = e = 0), which aUows for an exact 
solution of the problem. Indeed, the Hamiltonian has 
doubly degenerate eigenstates which can be represented 
by displaced oscillator states Q. In this degeneracy 
limit, one obtains two displaced harmonic oscillator wells, 
with equilibrium positions qo = ±X/muj^, so that the 
qubit can be localized in either the left or the right well. 
The system's wave function can be expanded in terms of 
a complete set of these orthonormal states and, for all su- 
perposition of the two degenerate levels, the ground state 
entanglement is found to increase with the coupling, be- 
coming maximal as A ^ oo. For A ^ (but e = 0), 
a numerical analysis is then performed in the same ba- 
sis, with some necessary truncation of the Hilbert space 
of the system. This approximation is known as "quasi- 
degenerate limit" and it has been shown to provide an 
accurate description of the system only for A/uj < 1, 
see Ref. JJj, due to the lack of orthogonality between 
different displacements. 

Here, we work in the opposite regime, and assume a 
fast qubit, A ^ w, to perform the adiabatic approxima- 
tion as described in the following section. Ground state 
entanglement is evaluated in section UTTl while section Hvl 
gives some concluding remarks. 



II. ADIABATIC APPROACH 

The Born-Oppenheimer approximation scheme can be 
followed more plainly by rewriting the Hamiltonian of 
Eq. (P) as follows 

H=^[Q^ + P''+ Da, + iW + LQ) a,] , (2) 

where the oscillator coordinates representation has been 
introduced, 

Q = ^(«^ + fl), P = i^{a^~a), (3) 

together with the dimensionless parameters D = 2A/cj, 
W = 2e/uj and L = 2\/Vrm^. 

The basic assumption of the well-known adiabatic ap- 
proximation is that the total wave function of a compos- 
ite system with one fast (the qubit) and one slowly (the 
oscillator) changing part can be written as: 



dQct>{Q)\Q)®\x{Q)) 



(4) 



The states \x{Q)) are the eigenstates of the "adiabatic" 
equation of the qubit part for each fixed value of the slow 
variable Q, 

[Da, + {W + LQ) a,] \xa{Q)) = E±(Q)\xa{Q)) , (5) 
which gives the eignvalues 

E±{Q) = ±E{Q) = ±^JD^ + {W + LQf . (6) 



The two eigenstates of Eq. (O can be written as 

\xim - -^(A-(Q)l+>-^+(Q)h» , (7) 

\xum - ^ {A+m+) + A-{Q)\~)) , (8) 

where |±) are the eigenstates of cr^ with eigenvalues ±1 
and 



W + LQ 



E{Q) 



(9) 



The subscripts I and u refers to the lower and to the upper 
effective adiabatic potentials felt by the slow oscillator, 
respectively. 



UuAQ) ^ ^ [Q^ ± EiQ)] 



(10) 



As we are primarily interested on ground state properties, 
we will concentrate on Ui from now on. 

A special case of interest is the one with the qubit 
working at degeneracy, W = 0. In this case one obtains a 
symmetric Hamiltonian with conservation of total parity 
(given by the joint transformation Q — > —Q and a^ — > 
—az)- Introducing the dimensionless parameter 



ill 
2D 



A2 



^A ' 



(11) 



one can show that for a < 1, the potential Ui^^'^{Q) 
can be viewed as a broadened harmonic potential well 
with its minimum at Q = and U^^'^{0) = —A. For 
a > 1 , on the other hand, the coupling of the oscillator 
with the qubit splits the oscillator potential producing a 
symmetric double well with the minima at 



= ±Qo = ±^Va2-l 



with 



ur'{±Qo) = -jU+i] . 



(12) 



(13) 



Qo is used as a kind of order parameter in Ref. j^l, in 
the limit D oo. 

For W 0, the symmetry is broken, and for this reason 
we refer to W as the asymmetry parameter. The form 
of lower potential for two different sets of parameters is 
shown in Fig.(^). 

Having obtained the state of the qubit, the last step 
in the adiabatic procedure is now to evaluate the ground 
state wave function for the oscillator, (jjoiQ): to be in- 
serted in Eq. to obtain the fundamental level of the 
coupled system. This wave function satisfies the one- 
dimensional time independent Schrodinger equation 



HadMQ) 



2 rfQ2 



UiiQ) ] MQ)^EoMQ), 

(14) 
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FIG. 1: The lower adiabatic potential for D = 10 and q = 2. 
The dashed line refers to the symmetric, W — 0, case (dashed 
line), while the solid line refers to VK = 1. The case of frozen 
qubit {W — D = 0) would have given a pair of independent 
parabolas instead of the adiabatic potentials Ui^u of Eq. HIO^ . 
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FIG. 2: Normalized ground state wave function for the oscil- 
lator in the lower adiabatic potential, for _D = 10 and a = 2 
and with W — (dashed line) and W — 0.1 (solid line). 



with (j)^{Q)dQ — 1, and where is the lowest eigen- 
values of the adiabatic Hamiltonian defined by the first 
equality. 

In Fig.©, the wave function (poiQ) is shown with D ~ 
10 and a = 2 for both the degenerate, W — 0, and a 
slightly asymmetric case, W = 0.1. Notice that even a 
very small value of the asymmetry parameter W gives 
rise to a wave function almost localized in the lower well. 

Given the function (j)a{Q), we can evaluate the reduced 
density matrix for the qubit and obtain the ground state 
entanglement. This is done in the next section. 



III. REDUCED QUBIT STATE AND 
ENTANGLEMENT 

The reduced density operator describing the qubit 
alone, when the overall system is in the ground state 



can be written as 

Po= / dQ\MQ)){MQ)\^ 7:il + bx<Jx + b,a^) 

(15) 

where b = (a) is the Bloch vector, whose non-zero com- 
ponents are explicitly given by the following integrals 

b^^-J_ ^oiQ)^jdQ, (16) 

and 

In Fig. (O and ^ , we show the dependence on the dimen- 
sionless quantity a — /2D of the ground state expec- 
tation values defined by Eqs. (|16|l and H17|) . respectively. 

It is easily seen that is different from zero only for 
W 0, while in the symmetric case the population is 
equally distributed between the states |±) of the qubit. 
This is due to the inversion symmetry of the adiabatic 
potential, which, for finite D, implies that the system 
does not localize in any of the wells and, consequently, 
that the state of qubit does not have a well defined value 

of (7^. 

In the case W = 0, the integrand of Eq. Hl()|) becomes 
the product of the squared ground state wave function 
and the square root of a Lorentz function centered at 
Q — 0. In the limit a ^ 0, this integral reduces to 
the normalization condition for the ground state wave 
function and thus b^ — —1- In fact, it is possible to 
show that, for small a, the main effect of the qubit is 
to renormalize the value of the oscillator frequency by a 
factor k — ^/l — a. As a result, the adiabatic ground 
state wave function for the oscillator is approximately 
given by 

</.o(Q)c. (^^y expj-^Q^I ^ (18) 

so that 

bx^-l + ^, for a « 1 . (19) 

For a 3> 1 the ground state wave function is located in 
spatial regions far from Q — and thus b^ — 0. To 
obtain an analytic estimation for large a, we can take as 
an approximate adiabatic ground state for the oscillator 
the symmetric superposition 

MQ)-^{MQ) + <P-iQ)} , (20) 

with 

cl^±{Q)^(-y exp{-^{QTQoA , (21) 
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FIG. 3: The dependence of the ground state expectation 
value hx ~ (ux) as a function of the parameter a, for D — 10. 
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FIG. 4: The 2-component of the Bloch vector as a function 
of a, for D = 10. 



where k' — 1 



is, again, a renormalization factor 



for the oscillator frequency. 
Taking the dominant contribution in Eq. 116() . one gets 



1 

a 



for a 3> 1 , 



(22) 



which we checked to be in very good agreement with the 
numerical solution. 

For W 0, bz becomes non-zero and decreases mono- 
tonically with increasing a with —1 as limiting value for 
a 3> 1. This is due to the fact that the contribution 
dominates in the qubit Hamiltonian in this regime, and 
therefore the qubit stays in the state |— ). 

Even in the presence of an asymmetry W, the x com- 
ponent of the Bloch vector continues to grow monotoni- 
cally from —1 to when a increases from zero, with only 
quantitative deviation from the W — behavior. 

We can, thus, summarize by saying that for small a, 
the state of the qubit is the lower eigenstate of , while 
for large enough a the qubit is found to be in the lower 
eigenstate of az- Between these two extreme cases, a 
cross-over occurs, which becomes a true, sharp transition 
for very large D and W — (see below). 

The knowledge of the qubit reduced density matrix al- 
lows us to evaluate the entanglement in the ground state. 
A quantitative measure of the entanglement between the 



qubit and the oscillator is given by the tangle [16| | , which, 
for globally pure states, is defined as 

T - 2 [1 - Tr(p2)] . (23) 

This is an entanglement monotone, giving t ~ for a sep- 
arable state and reaching t — 1 for maximally- entangled 
states. In our case, Eq. H23|l becomes 

T^l^bl-bl (24) 

The tangle is shown in Figs. ® and |2l for differ- 
ent values of the D and W parameters, as a function of 
the dimensionless quantity a. For any W ^ Q, the en- 
tanglement increases with increasing a before reaching 
a maximum value; after that, it decreases again to zero 
for a ^ 1. As stated above, this is due to the fact that 
the state of the system factorizes in this limit liW^Q. 
In the symmetric case, see Fig.((2j), the entanglement be- 
comes maximal as the coupling increases and, in the strict 
adiabatic limit D oo, the tangle becomes discontinu- 
ous at the critical value a = I and rapidly increases from 
zero to one when a > 1. 

This result has a simple analytic derivation that is eas- 
ily obtained from the thermal ground-state of the system. 
The reduced density operator describing the ground state 
of the qubit may be found by tracing out the oscillator 
variables from the thermal state 



Z{(3) ' 

^ oo. Here Zifi) 



and by taking the limit /3 —> 
is the partition function. The thermal density p pos- 
sesses the full symmetry as the Hamiltonian H and, if 
the ground state is non degenerate, the zero-temperature 
state coincides with the ground state of the system. It is 
important to stress that this is not generally true for de- 
generate ground states (as in the case A = e = 0). When 
a degeneracy arises, each individual ground state may 
not possess the same symmetries of the Hamiltonian. In- 
stead, they are always shared by the zero temperature 
state, which is just an equal mixture of all the possible 
ground states. This situation does not occur in our case. 

By rewriting the Hamiltonian |(2Jl as 

H = —p" + -—q^ + Aa, + (e + Xq) a, , (26) 
2m 2 



where 



q = 



(27) 



we see that the limit m ~* oo and muj'^ ~* const, is 
equivalent to neglect the kinetic energy of the oscillator. 
In this regime, thus, one gets 



Z{P) = Tr / dq{q\ 



\q) 



= 2 



dqe 2 9 cosh 



f3jA^ + {e + \qr 



(28) 
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In the first row, the states \q), employed to perform the 
trace, are just the positfon eigenstate of the oscillator. 

The thermal reduced density matrix of the qubit can be 
formally written in the form of Eq. (|15|) and is obtained by 
tracing out the thermal reduced density matrix over the 
oscillator degree of freedom. The temperature-dependent 
expectation values of <7x and ax, respectively, are, then 



bx = - 



b, = - 



, e + Xq 
dq . , , e 



sinh[/3A(g)], (29) 



A(g) 



sinh[/3A(g)], (30) 



where A{q) = + (e + Ai?)^. 

We focus, again, our discussion on the ground state 
(/3 oo). In this limit, the partition function may be 
evaluated by the steepest descent method. For e = 0, the 
integrand of Ea. H28|l is symmetric around q — 0. When 
/3 — > cx) and a < 1 this function has only a sharp max- 
imum at the origin, while for a > 1, the integrand has 
two sharp maxima at g = zLAy^a"^ — 1/A, symmetrically 
displaced around zero, where the function has a shallow 
minimum. In this limit, one easily obtain = and 



bx = 



-1, a<l; 
-1/a, a > 1. 



Then, for the tangle, one simply gets 



l/a2 



a < 1; 
a > 1. 



(31) 



(32) 



The first l/Z?-correction to this result can be obtained, 
for large a by using the expansion given in Eq. (|22|l to- 
gether with the definition of the tangle, Eq. (|24|) . These 
results are shown in Fig. 10, where the solid line is a 
plot of the tangle in the asymptotic regime (D — > oo). 

The same procedure used above can be carried out in 
the asymmetric case provided a value of q is found, such 
that 



A 



Xq 



v/A2 + (e + Ag)2 



(33) 



This equation has three nontrivial solutions. Within our 
saddle-point scheme, in the limit /3 — > oo, we must re- 
tain only the solution qm that corresponds to the lowest 
minimum of the potential. Therefore, we can write 



Xq„ 



A{q.m) ' 



A{q,n) 



(34) 



and, thus, one gets t — for any finite e 7^ 0. 

Indeed, it can be seen from Fig. that the tangle (for 
any value of a) decreases progressively with the increase 
of the asymmetry parameter W. Furthermore, as exem- 
plified in Fig. for any fixed non-zero value of W, the 
tangle approaches zero as D increases; so that, asymp- 
totically, the result implied by Eq. H34|l is obtained. 
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FIG. 5: The tangle r as a function of a for D = 10. Differ- 
ent curves, corresponding to diflFerent values of W, indicate 
the entanglement progressively decreases with increasing the 
asymmetry. 
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FIG. 6; The tangle r as a function of a for = 0.1 and 
different values of D. 



IV. CONCLUDING REMARKS 

In conclusion, we have discussed the adiabatic approx- 
imation for a qubit coupled to a single oscillator mode 
and we have derived the resulting entanglement in the 




FIG. 7: The tangle as a function of a in the symmetric case 
W = for different values of the qubit tunnelling amplitude 
D. One can appreciate that the result of Eq. (I32II is indeed 
reached asymptotically. 
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ground state, by giving simple analytical results in the 
strict adiabatic limit. The advantage of our approach, 
that requires a very small computational effort and cor- 
rectly describes the model system when A/w ^ 1, is 
to give a physically more transparent description of the 
ground state. 

As we have shown, the procedure is easily extended 
to the asymmetric case and this is important since the 
entanglement changes dramatically for any finite (how- 
ever small) value of the asymmetry in the qubit Hamil- 
tonian. As mentioned in section Ull above, this is due to 
the fact the this term modifies the symmetry properties 
of the Hamiltonian, so that the form of the ground state 
changes radically and the same occurs to the reduced 
qubit state. For example, for a large enough interaction 
strength, the qubit state is a complete mixture if = 0, 
while it becomes the lower eigenstate oi cr^ if W ^ 0. 
As a result, for large a, there is much entanglement if 
W = 0, while the state of the system is factorized and 
thus T = if ly ^ 0. This is seen explicitly in Fig. 
(O. Furthermore, from the comparison of Figs. l|njl, 
and Q, one can see that, with increasing a, the tangle 
increases monotonically in the symmetric case, while it 
reaches a maximum before going down to zero ii W 0. 

This is due to the fact that, in the first case, the ground 
state of the system becomes a Schrodinger cat-like entan- 
gled superposition, approximately given by 




for a > 1 , (35) 



where |^±) are the two coherent states for the oscillator 
defined in Eq. (|^ . centered in Q = ±Qo, respectively, 
and almost orthogonal if a ^ 1. 

In the presence of asymmetry, on the other hand, the 
oscillator localizes in one of the wells of its effective po- 
tential and this implies that, for large a, the ground state 
is given by just one of the two components superposed in 
Eq. (|35|) . This is, clearly, a factorized state and therefore 
one gets t — 0. 

Since r is zero for uncoupled sub-systems (i.e., for very 
small values of a), weather = or not, and since, for 
7^ 0, it has to decay to zero for large a , it follows that 
a maximum is present in between. 

In fact, for intermediate values of the coupling, there is 



a competition between the a-dependences of the two non 
zero components of the Bloch vector. In particular, the 
length l&l is approximately equal to one for both small 
and large a's, see Figs. ©-(El), but the vector points 
in the x direction for a <ti 1 and in the z direction for 
a 3> 1. The maximum of the tangle in the asymmetric 
case occurs near the point in which ~ 6^. 

For the symmetric case, we were also able to derive 
analytically the sharp increase of the entanglement at 
a = 1. This behavior appears to be reminiscent of the 
super-radiant transition in the many qubit Dicke model, 
which, in the adiabatic limit, shows exactly the same 
features described here, and which can be described along 
similar lines. 

Finally, we would like to comment on the relationship 
of this work with those of Refs. Q and The approach 
proposed by Levine and Muthukumar, Ref. employs 
an instanton description for the effective action. This 
has been applied to obtain the entropy of entanglement 
in the symmetric case, in the same critical limit described 
above. It turns out that this description is equivalent to 
a fourth order expansion of the lower adiabatic potential 
Ui. This approximation, although retaining all the dis- 
tinctive qualitative features discussed above, gives slight 
quantitative changes in the results. 

Concerning the asymmetric case, our results for the 
ground state entanglement appear similar to those found 
by Costi and McKenzie in Ref. , where the interaction 
of a qubit with an ohmic environment was numerically 
analyzed. It turns out that, for a bath with finite band- 
width, the entanglement displays a behavior analogous to 
that reported in Figs. (jSJ-JHI, when plotted with respect 
to the value of the impedance of the bath. Here, instead, 
we concentrated on the dependence of the tangle on the 
coupling strength between the qubit and the environ- 
mental oscillator. Unfortunately, the coupling strength 
is not easily related to the coefhcient of the spectral den- 
sity used in Ref. 0, and therefore one cannot make 
a precise comparison between the two results. At least 
qualitatively, however, we can say that the ground state 
quantum correlations induced by the coupling with an 
ohmic environment are already present when the qubit is 
coupled to a single oscillator mode. 
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